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The e expansion (expansion around four spacial dimensions) developed by Nishida and Son for 
a cold fermi gas with infinite scattering length is extended to finite scattering length to study the 
BEC-BCS crossover. A resummation of higher order logarithms and a suitable extension of fermion 
coupling in d-dimensions are developed in order to apply the theory in the BCS regime. The 
ratio between the chemical potential and the Fermi energy, fi/sp, is computed to next-to- leading 
order in the e expansion as a function of rj — l/(afcf), where a is the scattering length and kp 
is the Fermi momentum in a non-interacting system. Near the unitarity limit |?j| — > 0, we found 
m/ef = 0.475 - 0.707t; - 0.5t? 2 . Near the BEC limit 77 — ^ 00, fi/e F = 0.062/7? - if, while near the 
BCS limit n — > —00, (x/bf = 1 + O.7O7/77. Overall good agreement with Quantum Monte Carlo 
results is found. 



I. INTRODUCTION 

BEC-BCS crossover is a field attracting lots of attention recently 0, |E IE EL HI ■ The simplest systems to 
study the crossover are dilute Fermion systems with attractive interactions, for which the effective range 
for two-body scattering is much less than the inter-particle spacing. At zero temperature, such systems 
are characterized by a dimensionless number 77 = \/{akp), where a is the two-body scattering length, 
and kp is the corresponding Fermi momentum in non-interacting systems. For 77 large and negative 
(weak attraction) one finds the BCS solution with pairing and superfluidity. With 77 large and positive, 
corresponding to strong attraction with a two-body bound state well below threshold, the bound pairs 
will Bosc-Einstein condense (BEC). Experimentally. 77 can be chan ged between the BEC and BCS limits 
at will using the technique of Feshbach resonance [E IE IE IE HE HE U3 Theoretically physical quantities 
are expected to be a smooth function of 77, since there is no phase transition between the BEC and 
BCS limits. In each of these limits the behavior is nonperturbative, however, the effective interaction 
is weak; the system can be successfully described in a mean field approximation. On the other hand, 
dilute fermion systems near the unitarity limit (the 77 = limit) require treatments beyond the mean 
field approximation such as numerical simulations [IE HE HE HE HE HE LLLfl ■ 

Recently, Nishida and Son ha ve prop osed an analytical approach called "e expansion" to deal with 
problems near the unitarity limit |2[j |21| . In this approach one computes physical quantities in d spacial 
dimensions, then treats e = 4 — d as a small expansion parameter, and then sets e = 1 at the end to 
obtain the three dimensional physics. It is similar to the dimensional expansion method used by Wilson 
to study critical exponents of second order phase transitions |22j |. Nishida and Son's idea was inspired by 
the observation of Nussinov and Nussinov [23 that the ground state of a two-component fermion system 
in the unitarity limit is a free Bose gas for d > 4 (also see Refs. [3 HE f° r simplification of certain 
diagrams with d — > 00). It is because the two Fermion bound state (with zero binding energy) wave 
function behaves as $(r) cx l/r d-2 , where r is the separation between the two Fermions. The probability 
integral / d d r |<&(r)| has a singularity at r — > when d > 4. Thus the bound state has zero size and 
will not interact with each other. This feature was explicitly implemented in the theory set up near the 
unitarity limit in Refs. |2Ct l2ll|. As we will see, away from the unitarity limit, one also expects the 
localization of the bound state wave function all the way to the BEC limit but not to the BCS limit. 

Now the e expansion has been applied to few-body scattering [2(|, polarized fermions (i.e., with uneven 
chemical potential between two fermion species) in the unitarity limit p7^ . and near the unitarity limit 
[2ll| to identify the "splitting point" in the phase diagram where three different phases meet [23 • The 
critical temperature in the unitarity limit has also been computed |2^. Recently, the first next-to-next- 
to- leading order (NNLO) calculation in the e expansion was carried out for [i/sf in the unitarity limit, 
where fj, is the Fermion chemical potential and ep is the Fermi energy in non- interacting systems |30j |. 
Contrary to the nice convergence seen previously at next-to- leading order (NLO), the NNLO of fi/ep is 
as large as the leading order (LO) with an opposite sign. Even so, it is still premature to claim that 
the large NNLO correction cannot be tamed through reorganizing the series. Thus, we consider it is still 
worth while to explore the e expansion at lower orders. 

In this work, we apply the e expansion to NLO to study the BEC-BCS transition with an arbitrary 77. 
In the BEC side, an interesting relation between \x and the two fermion binding energy B is implemented: 
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2n + B = 0{e). (1) 

This comes from the chemical equilibrium between two fermions (with chemical potential 2[i) and one 
boson (with chemical potential —B) in a 4-d system where the bosons do not interact with each other. 
We also perform the resummation of the (e log r/) n jn\ terms arising from rf — 1. These logarithms are 
formally higher order in e but numerically important near the unitarity limit (e = 1, |rj| <C 1), and near 
the BEC and BCS limits (e = 1, |t;| 3> 1). For example, this resummation gives the linear term in (i/sf 
when <C 1: 

— = 0.475 - 0.707?7-0.5?7 2 . (2) 

Without this resummation, physical quantities are always even in r\ which would prevent us from reaching 
the BCS side. 



II. e EXPANSION 



We are interested in a two-component Fermi gas with zero range attractive interactions characterized 
by the scattering length a. In Nambu-Gor'kov formalism, the two-component fermions are denoted as 
"J" = ip\) T ■ After the Hubbard- Stratonovich transformation, the Lagrangian of the system is 

L = & (id t + — cr 3 + 11,73} * + tfa+Vd) + tfa-^d* - — , (3) 
V 2m J c 

where m is the fcrmion mass, Oi is the Pauli matrix acting on the spin space, and the coupling c depends 
on the two-body scattering length a. The ground state is a supcrfiuid state with a condensate (cf>) = 4> a , 
whose phase is chosen to be zero. <p is then expanded around <f>a, 

(f> = (f>o+ g<p- (4) 

The Lagrangian will be further decomposed into perturbative and non-perturbative parts. We will first 
look at the BEC side, than the BCS side. 



A. Power Counting in the BEC Side 

The Lagrangian can be decomposed as L = Lq + L\ + L-2, 
L = 



V 2 B 

«9 f +(—---) cr 3 + </> er+ + 



L 2 = -<p*(idt + ^— + —^\<p-fjt B ip*<p-g—(<p* + ip), (5) 
V 4m Co J c 

with B = 1/ (ma 2 ), /is = 2/i + B ~ 0(e) from Eq. which will be checked explicitly later. Near the 
BEC limit (77 > 1), B does not vanish in 4-d, thus B - O(e ) and fj, = [p B -B)/2 - O(e ). Similarly <j> 
does not vanish in 4-d in the BEC side. So <fio ~ O(e ). cq is determined by demanding the two fermion 
scattering amplitude has a l/(p — i/a) pole in momentum p in arbitrary dimensions |26( : 



c (47r) 2 - £ / 2 a 2 - £ V 2 

In 3-d (e = 1), 1/co = — m/(47ra) where a could be either positive (in the BEC side) or negative (in the 
BCS side). However in 4-d (e = 0), the scattering amplitude has a double pole l/(p 2 + 1/a 2 ). Thus, 
the sign of a does not play a role in 4-d and the difference between the BEC and BCS sides becomes 



ambiguous. Besides, if a < 0, 1/cq picks up a phase when < e < 1. To fix this phase, instead of using 
Eq. |(5J| , we will use 



— = signla] ■ -. 

co (47r) 2 - e / 2 |a| 2_e 

and write the \a\ dependence in 1/co as 

\a\- 2+e = a- 2 {l + eX), 



(7) 



X = (|<-l)/ e = £ 



(elog|o|) n 



(8) 



such that all the higher order logarithms are resummed and lumped into NLO. Then \a\ 2+6 — > 1/ \a\ at 
NLO after e — > 1 is taken at the end. Using this prescription, 1/cq ~ C(e _1 ). 

Lq is the non-perturbative part of the lagrangian which defines the fermion and boson propagators, 



Sf(p,B) 
D{p) 



-1 



{ Po (l + i0+)} z -E* 
-1 

p (l + J 0+)- £p /2' 

■ A2 



B/2 



Po-s p - B/2 



(9) 
(10) 



where e p = p /(2m), £^ = (e p + B/2) + 4>o and S/2 plays the role of effective mass for the fermions. 
L\ gives the perturbative interaction with g ~ ©(e 1 / 2 ) and /^s ~ ©(e 1 ), while L2 gives the counterterms 
for boson self energy and tadpole diagrams. The boson- fermion coupling g is tuned to 

6/4 
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•\/87r 2 e / m<^o 



2tt 



G( £ 1/2 ), 



(11) 



such that the boson self energy diagram without (with) an insertion of fis will be canceled by the 
first (second) term in Li in 4-d. The last term in L2 cancels the leading boson tadpole diagram when <po 
is chosen to minimize its effective potential. 

The quantity we will compute in this paper is p/ef- Following the procedure of Ref. [23], we first 
compute the effective potential V e s{4>o) and find the solution (f> that minimizes V c g. From V e s, the 
fermion number density is 



n f 



dV eS 



Op 



Then ep is 



The leading diagrams for the effective potential V s(<j>o) are listed in Fig. 1. 
The tree level diagram in Fig. 1(a) gives a contribution 



2tt 












TO 





(12) 



(13) 



B 



co WO J 



CO 



) = - = ^r -1 + - [ £ Q(0 O ) + 1] 



4>o m 
2ir 



2-e/2 



where 



sign[a\ 



72 



O(e ) 



(14) 



(15) 



The one- loop diagram in Fig. 1(b) yields 



1L [QO J 



«(2tt) 



Pom 
2tt 



5-£ 



lndet [S^faB)] 



*(*) 



l-e/2 



P / -2+e e . 9. -4<ftp 
2^1 1 4 1 4 j A 



22-e/2 S in(f )r(2-e/2) 



(16) 
(17) 



(a) 



(£>) 



(c) 




(d) 



FIG. 1: Feynman diagrams for the (f>o effective potential up to 0(e). The solid triangle is the 4>o/co potential. 
The solid lines are fermion propagators and dashed lines are boson propagators. The cross denotes one insertion 
of /i B . 



where 2-F1 is a hypergeometric function. The one-loop diagram in Fig. 1(c) with one insertion of \xb 
yields 



%b(w) = -Mb — ^ 0{1). 



Also, the two-loop diagram in Fig. 1(d) yields 

dp 5 ' e dq 5 - e 



i(27r) 5 - e i(27r) 5 - e 
4> m 



dB 



Tr [S F (p, B)a + S F ( qi B)aJ [ D(p - q) 



(18) 



-C(B/<M e^o 
0(e) , 



2tt 



(19) 



where C (2) is a decreasing function of z which converges numerically even at 4-d. The expression for 
C (z) and its numerical plot arc shown in the Appendix. So now we have 



VeS = V ca + V 1L + VlLjilB + ^2L- 



(20) 



B. Near the Unitarity Limit 

Now we study [ijep near the unitary limit by small B expansion of V e g to O(B). The 0(e~ 1 ) contri- 
butions cancel as 



TfP 1 dp 1 

-V C M = V 1 M = -^B + O(e ). 



The resulting 4> that minimizes V e s is 
0o = ^ + |(2Q(0 O ) + 1) 

^(3Ci + ln2 - 1) + B jc 2 + (3Ci + In 2 - ^ 



+ 



l\ 1 TT 1 



2 J 4 

where d = C (0) = 0.1442 [2fJ and C 2 = C (0) = -0.223 arise from V 2 l(</>o) with 



2tt 



•0(0- 



(21) 



e + 0(e 2 ),(22) 



(23) 



Since <fi ~ 0(e°), Eq. I|22[) shows that fiB should be 0(e), which agrees with the intuitive argument given 
above. 

Next, we obtain Ep as a function of \xb'- 

_ B ^ £l/2 _B_f c , Q&o) 



CSf Z£p 



48 "I 1 



0(e 5 / 2 ). 



(24) 



By taking the LO of <f> in Eq. (1221) and substituting [Ib by the LO expression in Eq. (|24(l . one obtains 
4>o = SFf 1 ^ 2 + 0(e) which implies sf ~ 0(e -1 / 2 ). </> can further be solved to higher orders in e by 
interaction. Finally, we obtain 



JL 

e F 



£ 3/2 

2 



elne 5(1 - In 2) - 120! 



5 



£^e 



1/2 



-e/2 



— +0(e 7 / 2 ,Be 5 / 2 ) 



0.475- 0.707?7~ 0.5r/ 2 , 



(25) 

where the arrow in the second line denotes extrapolation to three dimension e — > 1, and we have used 
B/ef — 2(akp)~ 2 = 2rj 2 . Note that the result of Ref. [2l| is reproduced if we truncate the expansion 

of (B/e F e 1/2 y as 1 - e/21og[B/e F ] + e/41og[e], which does not give a linear term in r\. In order to 
compare with the Quantum Monte Carlo (QMC) results [l3l IT^j. /i/ef can be converted to energy per 
particle relative to the value of a non-interacting fermi gas 

E/A 



E /A 



(26) 



with E /A = 3/5e F in 3-d. The relation is, fx/e F = £-f f^- Our result gives £ ~ 0.475-0.884r;-0.833r ? 2 . 
This is compared with two sets of QMC results denoted as QMCl[l3| and QMC2 Q in Fig. 2. 



The BEC Limit 



In BEC limit, B » ef or 77 > 1, we perform 1/B expansion to 0{B~ 2 ). Note that V2l(<M ~ 0(B~ 2 ) 



and is neglected. The effective potential has a simple expression 



V'c 



eft [<P0 



B(Bmy/ 2 sin(f) T (1 



0(B~ 2 ), 



(27) 



which is minimized when 4> = <f> Q = 2^/2B[LbI^- Following the same procedure, we obtain 



JL 

e F 



B 
~2e~F~ 
O.O6277- 1 - if 



5< 



3 -la 4 



1/2-1 



(28) 



This result is corresponding to £ = 0.052?7 1 — |?? 2 . Note that the rj 2 term just comes from the free 
molecule contributions, while the vj~ 1 correction is from the §\ term in V e ff, corresponding to the mean 
boson-boson interaction. The numerical value agrees well with the QMC1 result, which yields £ = 
O.t^rr 1 - f?? 2 0. The QMC2 calculation used a finite range cosh shape potential, thus the deviation 
from the QMC1 curve indicate the finite range effects. 

D. Power Counting in the BCS Side 



Next we will turn to the BCS limit, —77 3> 1, where the situation changes drastically. As —a — > 0, 
the gap vanishes and the ground state should be a non- interacting Fermion gas with A = 4> /fJ, <C 1. 
Thus we will count /1 to be 0(e ). And since Eq. (JJJ is no longer valid, we will write the Lagrangian as 
L = Lq + L\ + L2, where 



L 
Li 
Li 



id t + 
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t> a + 



t> <7- 



* + ip* [id t + — +2n 
4m 



Co 



idt + I^ + 2 ^) 9~^* + 



(29) 
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FIG. 2: (Color online) £ (energy per particle relative to the value of a non-interacting fermi gas) shown as 
functions of (akp)^ 1 . QMC1 flal and QMC2 [T3[ are Quantum Monte Carlo calculations with zero and finite 
range, respectively. NUL is our near unitarity limit result of Eq. 1251 in the e expansion. BEC is our result in 
Eq. 1281 . BCS is obtained from numerical solving Eq. 13611 . The lower panel is the blow-up of the near unitarity 
regime. 



In Lq the boson has an "effective mass" g 2 /co = B + 0(e) which will suppress the boson propagation 
in the BCS limit where B = 1/ma 2 is big. Again, the boson self-energy renormalization vanishes at 
4-d. The effective potential receives leading contributions from Fig. 1(a), (b), and (d) and yields 

lndet [S/(p,-2/i)] 
l{2 ^- H ^ f -J T [^fo ~^WS F {q, -2/x)cr_] D(p - q), (30) 

where Sf(j>, — 2/i) is the one defined in Eq.© with the substitution B — > — 2/i. The boson propagator 
D^ 1 {p) = — J?o(l + iS) + Sp/2 + B — 2/x has been replaced by £) -1 (p) = B which is a good approximation 



when B 3> 2fi or cquivalcntly — 77 ^> 1. Since we expect A <C 1, by expanding to C(A 2 ) we obtain 



1 



12tt 2 



l + £ (3^/2-1 +ln ( 2 ^)_ 7 + 



11 



8tt 2 



gl-e/2 _ 2 



1 - 7 + In (47r) - ! f/9 ~ /9 , , A 2 
^ — LB 1 -*' 2 + 2B e/2 ~ 1 + 7 + In 



8tt 



(31) 



where -B = -B//i. 

At this order, the gap equation, dV c s/dA = 0, has a simple solution 



elnA = 1 e 



ln(47r) zl±± B | g- 1 



1 + 7 - ln(87r) 



0(e 2 ,A 2 



(32) 



Thus, In A ~ C(e 1 ). We will formally count A ~ 0(e 1 ^ 2 ) in Eq.JjHJ to indicate the smallness of the 
gap. 

The rif and ep can be computed using Eqs. 1)121) and 1)13)1 : 



JL 



1 - eB 6 / 2 " 1 



3A 2 



B 



2e 

67-6 ln(87r) 
72~ 



1 



£l-e/2 



11 



1 - 2 7 + 10 In 2 + 2hi7r 



In A 



Since B 



2 £F 
\? 



(akp) M 



and fx/sp = 1 + £>(e°), Eq. (gU shows that 



A ~ exp 



1 — \aki 



l-2+e 



(33) 



(34) 



This implies that, in 4-d, the gap vanishes and y^jep = 1 when —akp <C 1. In 3-d and with —akp <C 1, 
the gap is exponentially suppressed. One might argue that A should be smaller than 0(e 1 ^ 2 ) as counted 
in Eg. 1)31)) . This is certainly true when e — > 0. However, since we will take e — > 1 at the end, this counting 
suits our purpose especially when working between the BCS and unitarity limits. 

To compare with the low energy theorem £ = 1 + + C(?7~ 2 ) obtained in 0,112, we perform a 

I/77 expansion to Eg. 1)33)) . At NLO (0(e)) with all the exponentially suppressed C(A 2 ) terms neglected, 
we see no rf 1 correction. However, we expect that at NNLO, the e.B e / 2-1 will be fully recovered and 
gives 



£p 



akp 

7f 



(35) 



as e — > 1. This implies £ = 1 + 0.589^ _1 + 0{rf~ 2 , e 2 ). This agrees reasonably well with the low energy 
theorem. 

As — r) becomes smaller but not much smaller than unity, the C(A 2 ) effect in [i/sf becomes more 
important. From the thermodynamic relation E/A = V c s/nf + fj,, we obtain 



27 B~ l + [37 - 3 ln(27r) - 4] B - 13 - 67 + 6 ln(27r) 



A . 



(36) 



Eqs. 1)32)) and 1)36)) give £ in terms of B, while Ea. ()33)l gives B in terms of akp. The numerical result of 
£ as a function of akp is shown as the BCS curve in Fig. 2. 



III. CONCLUSIONS 



We have extended the e expansion (expansion around four spacial dimensions) developed by Nishida 
and Son for a cold fcrmi gas with infinite scattering length to finite scattering length to study the BEC- 
BCS crossover. A resummation of higher order logarithms and a suitable extension of fermion coupling 
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FIG. 3: The C(z) function that appears in Eq. H19H 

in d-dimcnsions (see the discussion around Eqs. GEO have been developed in order to apply the theory 
in the BCS regime. The ratio between the chemical potential and the Fermi energy, h/ef, has been 
computed to next-to- leading order in the e expansion as a function of r\ — \/{akp). Near the unitarity 
limit \rj\ -> 0, wc found fi/e F = 0.475 -0.7077j-0.5r; 2 . Near the BEC limit t] -> oo, fx/e F = 0.062/r; - if, 
and near the BCS limit r\ — > — oo, /i/ep = 1 + 0.707/ry. As shown in Fig. 2, overall good agreement with 
Quantum Monte Carlo results has been found. 




IV. APPENDIX 

The function C (z) that appears in Eq. I|19|l is defined as 

(Vx 2 + 1 - x) (Vy 2 + 1 - y) 



C{z) 



oo p oo 



dx I dy- 

z/2 Jz/2 



2Vx 2 + ly/y 2 + 1 



-z + x + y + 2\Jx 2 + 1 + 2^y 2 +1 - J -z + x + y + 2 (y ' x 2 + 1 + s/y 2 + lj - (z - 2x) (z - 2y) 



The numerical value of C (z) is plotted in Fig. 3. 

In Ea. l3Ufl . the one fermion loop contribution to the effective potential near the BCS limit has the 
analytic expression: 



dp 



5-e 



i{2n) 



AT 1- 



5-e 



lndet [^- 1 (p,-2 M )] 
e-6 



21- 



^A 2 



2F1 



r(2 

e-6 e-2 1 

4 ' 4 '2' 



1 

A 2 



2r 



- 1 



2-fl 



e e 
4' 4 



2' A 2 
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